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Abstract. Trace formulae for d-regular graphs are derived and used to express the 
spectral density in terms of the periodic walks on the graphs under consideration. The 
trace formulae depend on a parameter w which can be tuned continuously to assign 
different weights to different periodic orbit contributions. At the special value w = 1, 
the only periodic orbits which contribute are the non back- scattering orbits, and the 
smooth part in the trace formula coincides with the Kesten-McKay expression. As 
w deviates from unity, non vanishing weights are assigned to the periodic walks with 
back-scatter, and the smooth part is modified in a consistent way. The trace formulae 
presented here are the tools to be used in the second paper in this sequence, for showing 
the connection between the spectral properties of ci-regular graphs and the theory of 
random matrices. 



1. Introduction and preliminaries 

Discrete graphs stand at the confluence of several research directions in physics, 
mathematics and computer science. Notable physical application are e.g., the tight- 
binding models which are used to investigate transport and spectral properties 
of mesoscopic systems [1], and numerous applications in statistical physics (e.g., 
percolation [2]). The mathematical literature is abundant with studies of spectral, 
probability and number theory, with relation to discrete graphs [32, 3, 33]. Models of 
communication networks or the theory of error correcting codes in computer science use 
graph theory as a prime tool. 

In the present series of papers, we would like to add yet another link to the list above, 
namely, to study graphs as a paradigm for Quantum Chaos. In making this contact we 
hope to enrich quantum chaos by the enormous amount of knowledge accumulated in 
the study of graphs, and offer the language of Quantum Chaos as a useful tool in graph 
research. 

The first hints of a possible connection between Quantum Chaos and discrete graphs 
emerged a few years ago when Jakobson et.al. [4] studied numerically the spectral 
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fluctuations for simple (i-regular graphs (graphs where the number of neighbors of each 
vertex is d and no parallel or self connections are allowed). In particular, they sampled 
randomly the ensemble of rf-regular graphs, computed the spectra of the adjacency 
matrices, and deduced the mean nearest-neighbor distributions (for d = 3,4,5). They 
found that within the statistical uncertainty, the computed distributions match the 
prediction of Random Matrix Theory. The work of Terras [34] should also be consulted 
in this context. The purpose of the present work is to adopt the techniques developed in 
Quantum Chaos to investigate the connection between the spectral statistics of d- regular 
graphs and Random Matrix Theory. The main tools which we shall use to this end are 
trace formulae, which, in the present case, relate spectral statistics to the counting 
statistics of periodic walks on the graphs. In the present paper, the first in this series, 
we shall develop the tool kit - namely - will derive trace formulae for regular graphs. We 
shall show that a large variety of trace formulae exist, all of them provide expression for 
the same spectral density but using differently weighted periodic walks. We shall also 
show that there exists an optimal trace formula, in the sense that the smooth density 
coincides with the mean density (with respect to the Qv,d ensemble), for this optimal 
trace formula, the oscillatory part stems only from a subset of periodic orbits. These 
are periodic orbits in which there are no back-scattering (reflections). In the second 
paper in the series, wc shall use the optimal trace formula to obtain some results which 
support the conjecture that spectral statistics for regular graphs follow the predictions of 
random matrix theory. In many respects, the present study follows the development of 
the research in Quantum Chaos where it was conjectured [5] that the quantum spectra 
of systems whose classical analogues are chaotic, behave statistically as predicted by 
Random Matrix Theory. This conjecture, which was originally based on a few numerical 
studies, brought about a surge of research, and using the relevant (semi-classical) trace 
formula [6], the connection with Random Matrix Theory was theoretically established 
[7, 8, 9]. 

To provide a proper background for the ensuing discussion, we have to start with 
a short section of definitions and a summary of known facts. 

1.1. Definitions 

A graph ^ is a set V of vertices connected by a set S of edges. The number of vertices 
is denoted by V = \V\ and the number of edges is E = \S\. The V x V adjacency 
[connectivity) matrix A is defined such that Aij = s if the vertices i,j are connected by 
s edges. In particular, j = 2s if there are s loops connecting the vertex i to itself. A 
graph in which there are loops or parallel edges, is called a multigraph. 

In the present work we mainly deal with connected simple graphs where there are 
no parallel edges (Ajj e {0,1}) or loops [Ai^i — 0). The degree di {valency) is the 
number of edges emanating from the vertex, d^ — 'Ylij=i ^ ^-regular graph satisfies 
di = d i : 1 < i < V, and for such graphs dV must be even. The ensemble of all 
d-regular graphs with V vertices will be denoted by Gv,d- Averaging over this ensemble 
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will be carried out with uniform probability and will be denoted by (■■■). 

To any edge b = one can assign an arbitrary direction, resulting in two directed 
edges, e = [i.j] and e = [j, i]. Thus, the graph can be viewed as V vertices connected by 
edges b = 1, ■ ■ ■ , E or hj 2E directed edges e = 1, ■ ■ ■ , 2E (The notation b for edges and 
e for directed edges will be kept throughout). It is convenient to associate with each 
directed edge e = [j, i] its origin o(e) = i and terminus t{e) = j so that e points from 
the vertex i to the vertex j. The edge e' follows e if t{e) = o(e'). 

A walk of length t from the vertex x to the vertex y on the graph is a sequence of 
successively connected vertices x — Vi,V2, • ■ ■ ,Vt — y. Alternatively, it is a sequence of 
t — 1 directed edges ei, ■ ■ ■ , et-i with o{ei) — Vi, t{ei) = Wj+i, o(ei) = x, t(e^_i) — y. A 
closed walk is a walk with x = y. The number of walks of length t between x and y 
equals (A*)^ .j,. The graph is connected if for any pair of vertices there exists a t such 
that {A%^, ^ 0. 

We have to distinguish between several kinds of walks. There seems to be no 
universal nomenclature, and we shall consistently use the following: 

A walk where Cj+i T^e^, , 1 < i < t — 2 will be called a walk with no back-scatter 
or a nb-walk for short. 

A walk without repeated indices will be called a path. Clearly, a path is a non 
self-intersecting nb-walk. 

A t-periodic walk is a closed walk with t vertices (and t edges). Any cyclic shift 
of the vertices on the walk produces another t-periodic walk (which is not necessarily 
different from the original one). All the t-periodic walks which are identical up to a 
cyclic shift form a t-periodic orbit. A primitive periodic orbit is an orbit which cannot 
be written as a repetition of a shorter periodic orbit. 

Amongst the t-periodic orbits we shall distinguish those which do not have back- 
scattered edges and refer to them as periodic nb-orbits. The frequently used term cycles, 
stands for periodic paths (non self intersecting nb-orbits). 

In order to count periodic walks, it is convenient to introduce the 2E x 2E matrix 
B which describes the connectivity of the graph in terms of its directed edges: 

Be,e' = St(e),o{e') ■ (l-l) 

The matrix which singles out edges connected by back-scatter is 

Je,e' = Se,e' ■ (1-2) 

The Hashimoto connectivity matrix [10] 

Y^B-J, (1.3) 

enables us to express the number of t periodic nb- walks as try*. A slightly more general 
form: 

Y{w) ^ B - wJ , weC (1.4) 

gives a weight 1 to transmission and weight 1 — w to back-scatter. Now, try*(w) = 
'^„N{t; g){l — wy, where N(t;g) is the number of t periodic walks with exactly g 
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back-scatters. Clearly, tr Y*{w) can be considered as a generating function for counting 
periodic walks with specific t and g: 



N{t;g) 



g\ dw^ 



(1.5) 

w=l 

The discrete Laplacian on a graph is defined in general as 

L=-A + D, (1.6) 

where A is the connectivity matrix, and D is a diagonal matrix with Dj^j = di. It is 
a self-adjoint operator whose spectrum consists of V non negative real numbers. For 
d-regular graphs D is proportional to the unit matrix and therefore it is sufficient to 
study the spectrum of the adjacency matrix A. This will be the subject of the present 
paper. 

The spectrum a {A) is determined as the zeros of the secular function (characteristic 
polynomial) 

Z^(/x) = det(/x/(^) - A) . (1.7) 

Here, is the spectral parameter and I^^^ is the unit matrix in V dimensions. The 
largest eigenvalue is d, and it is simple if and only if the graph is connected. If the 
graph is bipartite, —d is also in the spectrum. 

The spectral measure (spectral density) is defined as 

Pit') = ^ E ^if^-f'a). (1.8) 

The "magnetic" Laplacian [11] is defined by 

Ljjf ) = - A,,e^^-- + rfAj , = (1-9) 

The phases attached to the edges (i,j) play the role of "magnetic fluxes" (to be 
precise, the phases are the analogue of gauge fields, and a sum of gauge fields over a cycle 
is a magnetic flux). The Laplacian is complex hermitian, and therefore the evolution it 
induces breaks time reversal symmetry. Again, for d-regular graphs it suffices to study 
the magnetic adjacency matrix 

= A.e^^^'^- . (1.10) 

The ensemble of random magnetic graphs consists of the random graph ensemble Qv,d 
with independently and uniformly distributed magnetic phases (pij. 
A weighted Laplacian is defined by: 

L^^^ = -AjW^j + W,,A,, , Wij = Wj,i WijeR. (1.11) 

Where Wij are weights defined on the edges of the graph. We shall restrict our attention 
to the weighted adjacency matrix: AfJ'^ = AijWij . The ensemble of random weighted 
graphs consists of the random graph ensemble Qv,d with independently and uniformly 
distributed weights Wij in the interval \Wij\ < 1. 
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1.2. Background 

Adjacency matrices of random d-regular graphs, have some remarkable spectral 
properties. An important discovery which marked the starting point of the study of 
spectral statistics for d-regular graphs, was the derivation of the mean spectral density 
by Kesten [12] and McKay [13]: 

/>™(/<) = Jim {p{p)) = { . (1.12) 

t for \n\ > 2^/d^ 

The proof of this result relies on the very important property of random d-regular 
graphs, namely, that almost surely every subgraph of diameter less than \og^_i K is a 
tree. Counting periodic orbits on the tree can be done explicitly, and using the close 
relations between these numbers and the spectrum, one obtains (1.12). 
If the entire spectrum of A (except from the largest eigenvalue) lies within the support 
[—2\/d— l,2\/d— 1], the graph is called Ramanujan (For a review, see e.g., [23] and 
references cited therein). 

Trace formulae for regular graphs were discussed in the literature in various 
contexts. The late Robert Brooks [14] studied the connection between the number 
of closed paths on a graph and its spectrum, and proposed a Selberg-like trace formula. 
Some of the ideas developed in the present work are related to Brook's results. In this 
paper we derive two families of trace formulae which depend on a continuous (complex) 
parameter w. The optimal trace formula alluded to in the introduction is obtained 
when w = 1. In this case the trace formula uses the subset of nb-walks which have 
some advantages (see for example [15, 16]). We shall see that with this choice of w, we 
get the Kesten-McKay measure as a limiting distribution, and an oscillatory part which 
vanishes under an ensemble average. We will make use of this part in the following 
paper in this series of papers. The trace formula for this choice of w coincides with a 
trace formula derived by P. Mnev, ([41]) in an entirely different way. 

A different, more physical approach, in which a scattering formalism was used, 
resulted in a trace formula [17] (see also [24]), which is formally similar to the trace 
formula for the spectrum of the Laplacian on metric graphs [19, 18]. For d-regular 
graph it reads 

\ di 1 ci ^ 1 

P(l^) = - 9 , ^9 + 7^ Jm — V - tr ([/(// + ie)*) (1.13) 



t=l 



where t/(Ai), the graph evolution operator, is a 2E x 2E unitary matrix defined as 



d — ifi 



B + J 



(1.14) 



d — ifj, d — ifi' 

Y, J and B were introduced before. The infinite sum in (1.13) can be written as a sum 
over ^-periodic walks. Denoting by g the number of back-scattering along the walk, 

trC/*= ^e^*(=^'-^*^'^^^-5)yAr(i;o) ((^ " 2)^ + /^^) ^ ^-.^ arctan ^ (^ -^^. 

(^2 + ^2)1 a 25 
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The first term in the trace formula, (as usual, an algebraic function of the spectral 
parameter), is referred to as the "smooth" spectral density. It consists here of a 
Lorentzian of width d centered at 0. The sum over the periodic walks is formal and 
it diverges at the spectrum of the adjacency matrix. We shall show below that there 
exists another continuum of trace formula and that (1.13) is obtained as a special case. 

It was already observed in [17] that the trace formula above is not satisfactory 
in the sense that the leading Lorentzian is very different from the asymptotic mean 
spectral density (1.12). It was suggested that a re-summation of the infinite sum 
would extract the difference between the Lorentzian and the Kesten-McKay expression. 
Indeed, summing up the contribution of the shortest, t=2 periodic orbits, one obtained 
a correction term which exactly canceled the ^ tails of the Lorentzian. However, a 
systematic re-summation of (1.13) to extract pMKif^') requires a new approach which 
will be presented in section (2). The discussion of this problem leads naturally to a more 
general question: Can one distinguish amongst the periodic orbits on the graph distinct 
subsets, each responsible to a different feature in the spectral density? The Kesten- 
McKay theory seems to favor an affirmative answer, since the asymptotic density is 
derived from tree hke periodic orbits. But what about the rest? Does one really need 
all the legitimate periodic orbits on the graph, or can one do with a subset? Is this 
distinction unique? These, and other questions, will be further discussed in the following 
sections. 

The paper is organized as follows. The Bartholdi identity is the corner-stone of the 
theory presented in the present paper. However, since the derivation is strictly technical, 
we leave it to an appendix. The appendix reviews shortly the Bartholdi identity and 
its proof, which is generalized here to include the cases of magnetic adjacency matrices, 
multigraphs and weighted graphs as well. 

In the first section, we derive the family of trace formulae, which depend on real w. 
Then, we specialize to the choice w = 1 both in the magnetic and non-magnetic cases. 
The different forms of the trace formula will help to elucidate the question posed in the 
previous paragraph regarding the different roles played by periodic orbits of different 
topologies. Furthermore, the dependence on w can be exploited to unravel combinatorial 
information about the graph at hand. Finally, the trace formula (1.13) will be shown to 
be but one of a (continuous) family of trace formulae which have one feature in common, 
namely that they are based on traces of unitary operators of the type (1.14). This is 
achieved by allowing w to assume complex values. 

2. A continuous family of Trace Formulae 

In the present section we shall derive a continuous family of trace formulae which depend 
on a real parameter w. The parameter w controls the weights which are given to periodic 
walks with different numbers of back-scatters. Consider the matrix Y{w) (1.4). For d- 
regular graphs, the Bartholdi identity reads (see Appendix A): 

det(/(2^) -s{B-wJ)) = {l- w^sY'^ det{I^^'\l + w{d - w)s^) - sA) . (2.1) 
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Its importance in the present context comes from the fact that it connects the spectrum 
of the adjacency matrix A with that of the matrices Y{w) — B — wJ, which can be used 
to count various types of cycles and walks on the corresponding graph. It implies that 
the spectrum of Y{w) = B — wJ is 

a{Y{w)) ^{{d-w),w, +w X {E -V),-w X {E -V), 

{^Jw{d-w) e^*^ ^Jw{d - u)) k^l,---{V - 1))} (2.2) 



where 0^ = arccos 



for all k 



2y^w{d — w) 

The fikS with k = 1, - ■ ■ {V — 1) are the non trivial eigenvalues of the adjacency matrix, 
whose spectrum is ordered as a non increasing sequence 

d — Ho > Hi > IjL2 > ■ ■ ■ > Hv-i > —d . (2.3) 

{Hv-i assumes the value —d if and only if the graph is bipartite. The bipartite graphs are 
rare in Qv,d and are excluded from the discussion from now on). If w is in the interval 
[1, ^^], we have 2\/d — 1 < 2^Jw{d — w) < \/d^ — \. This ensures that for generic 
graphs one can always find a value of w G [1, so that for all /c, |/ifc| < 2^w{d — w), 
and all the ^^'s are real. The freedom to choose w allows us to use the trace formula 
for almost all d-regular graph, and in particular for non-Ramanujan graphs. 
It is convenient to introduce the quantities yt{w), 



yt{w) 



1 tTYHw) 



{d 



w 



V {^w{d-w)y 

which (unlike trF*(w)) are bounded as t 
eigenvalues of Y{w) are used now to write, 

1 



(2.4) 



oo. The explicit expressions for the 



yt(w) = 



w 



V \d-w 



2 d-2 
+ — — 



w 



d — w 



v-i 



Hk 



y 2\iw\d — w) 



(2.5) 



k=l 



where Tt{x) = cos (t arccos x) are the Chebyshev polynomials of the first kind of order 



t. The fact that the yt{w) are bounded, is guaranteed since 



d—w 



< 1 whenever 



w e [1, ^^], and d > 3, and since the Chebyshev polynomials are bounded. 



Multiplying both sides of (2.5) by ::^^T,{^-^J^^: 

get: 

1 °° 1 



and summing over t, we 



t=o 



w{d 



w) 



oo 



U) 



(d - 2) 1 

27r -^(l + M 



+ 



1 



k=l 



2^w{d 




l^ 



'2-sJw{d — w) 



) + 



(2.6) 
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Where Srix^y) is defined by: 

2 °° 1 

St{x, y)^-Yl T^WT,{y) . (2.7) 

ST{x,y) is the unit operator in the L'^[—l, 1] space where the scaler product is defined 
with a weight :y5=f- Indeed, 

f_^^^^8T{x,y)f{x)^f{y). (2.8) 
For t = 0, 1, 2 one can easily show that 



, 1 -I d-w dil-wf Id-w 

yo = d- — ; yi = —\ ; 1/2 = —j-, r - ■ 

V V \l w w[d — w) V w 

Writing \/l — • S{x — y) = 5T{x,y), and using elementary identities involving the 
Chebyshev polynomials, we get an expression for the density of states which is supported 
in the interval |//| < 2-\Jw{d — w) : 

p^l,) = p^^ooth^,,; w) + p''-(/.; w) + w) , (2.9) 

where 



psmooth^^. ^) 



d/(27r) 



■\/Aw{d — w) — iJ? 



^ _ {d-2w){d-2) ^ {w-lf{^Ji^-2w{d-w)) \ 
— jj? w^{d — wY J 



1 / v-^ Vtiw) itarccos 

exp 



.t=3 



^/^w{d^^vo)^^~^ 



2y/w(d-w) 



P-M - , fl + ^ + l^-2wid-.) ^ d_2u,\ 

2'K^Aw{d-w) - ^J? V W W2 0?-/xy 

Equations (2.9,2.10) are the main result of this section. In form they are very similar 

to well known trace formulae from other branches of Mathematical Physics. They are 
composed of a smooth part p«"»°o*'^(^; u;) which is an algebraic expression in as shown 
in Figure (1), and an oscillatory part, computed from information about the periodic 
walks. The amplitudes are obtained from the (properly regularized) count of t-periodic 
walks, and the phase factors explicitly given in the second line, are the analogues of the 
"classical actions" accumulated along the walks. (2.9) is a generalization of (1.13), as 
will be shown in (4) 

It is important to notice that the left hand side of (2.9) does not depend on w, 
because the density of states of the adjacency operator depends only on the graph and 
not on the choice of w. Therefore, the right hand side must also be w-independent. In 
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0.04 - I ,■ I 1:1- 

^jj , , , , , ^ ^ W .1 

-6-4-3-2-1012345 

Figure 1. w) for 5 regular graphs for various values of w. 

solid: w = 1, dashed: w — 1.2, dotted: w — 1.5, dash-dot: w — 1.7. 



other words, the w dependence of the smooth part, is offset by a partial sum of the 
oscillatory part. This is reminiscent of the partition of the trace formula (1.13), where 
the smooth part is a Lorentzian, and it was shown that the sum over 2-periodic orbits 
gives a contribution which exactly cancels the leading 1/yU^ behavior for large 

Having the freedom to choose w, it is natural to look for the most appropriate 
or convenient partition of the spectral density into smooth and oscillatory parts. We 
shall show in the sequel that this is obtained when w = 1, because the smooth density 
coincides with the mean density (with respect to the Qv,d ensemble) (1.12). 

Finally, we also mention that a trace formulae for multigraphs can be derived in an 
analogous fashion. The smooth part remains unaltered, and in the oscillatory part, the 
sum starts from t = 1 since loops are allowed. 

2.1. Trace formula in terms of periodic nb-walks (w = 1) 

The case w = 1 plays a special role in the present theory and its applications. The fact 
that the smooth part of the trace formula is identical to the Kesten-McKay expression 
was mentioned above, and will be discussed further in the sequel. As will be shown, this 
is a direct consequence of the fact that the counting statistics of t-periodic nb-walks in 
the Gv,d is Poissonian for t < \og^_i V with (tr F*) = {d — 1)* [31]. 

The trace formula can be obtained by substituting w = 1 m (2.10). Alternatively, 
one can start from the Bass formula for d-regular graphs [21]: 

det(/(2^) - sY) = (1 - s^f-^ det(j(^)(l + (rf - l)^') - sA) . (2.11) 

where Y is the Hashimoto matrix defined in section (1), and trF* counts the number 
of t-periodic nb-walks. One can now follow the same steps as in the previous section. 
However, this requires restricting the discussion to Ramanujan graphs only. Under this 
condition, < 2\/ d — 1, k = 1, - ■ ■ {V — 1), and the spectrum of Y is 

a{Y) = {[d-l),l, +lx{E-V), -lx{E~V), 
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(Vd^ e''t'\ y/d^ e-''t"', A; = 1, • • • - 1))} (2.12) 

where 0^ = arccos — , < 0^ < vr. 
2\Jd — 1 

For large t, the number of t-pcriodic walks is dominated by the largest eigenvalue, so 
that asymptotically trF* ~ (d — 1)*. Shortening the notation by using yt — yt{w — 1) 
we have 

which is the properly regularized number of t-periodic nb-walks. Going through exactly 
the same steps as in section (2), we get: 

The first term is the smooth part, and can be identified as the Kesten-McKay density. 
For large t, since trF* is dominated by {d — 1)*, it is clear that in the limit of large 
V, yt tends to zero. The counting statistics of t-periodic nb-walks with t < log^_i V 
is Poissonian, with (trF*) — (d — 1)*. Thus, the mean value of yt vanishes as O [y). 
Hence, 

lim (p(/x)) = pxm(/^). (2.15) 

The above can be considered as an independent proof of the Kcstcn-McKay formula 
(1.12). The original derivation relied on the fact that regular graphs look locally like 
trees, for which the spectral density is of the form (1.12). Here, it emerged without 
directly invoking the tree approximation, rather, it appears as a result of an algebraic 
manipulation. 

The last term, ^ p'^°^'^ (fj,) , is the correction to the smooth part, for finite V, which 
results from yo,yi,y2- It is given explicitly by 

^ WWk^^ (l . . ^ - 2(. - 1) + 1^) (2.1C, 

Equation (2.14) is identical to a trace formula derived by P. Mnev ([41]). Our derivation, 
however, follows an entirely different path which allows us to place the trace formula 
as a special case in a more general setting, and provides the leading order correction, 
fp--(/x). 

If the graph at hand is non-Ramanujan, the theory should be modified since some 
of the phases (pk in (2.12) become complex. Consequently, the yt diverge exponentially 
with t and the resulting trace formula (2.14) is ill defined. There are two ways to 
circumvent this problem. Either to use the trace formula with > 1, at the cost of 
using periodic walks with back-scatter, or to remain with nb-walks but reformulate the 
trace formula so that it describes a coarse-grained version of the spectral density. We 
shall discuss below the latter option. 

Consider a non-Ramanujan graph Q and denote by 7(^) the set of eigenvalues 
which lie outside the Kesten-McKay support [—2\/d — l,2\/d — 1]. The rigorous and 
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numerical information available to date on the set of eigenvalues jiG), suggests that 
generic {V, d) regular graphs have the following properties [42, 43]. 

• The distance between ^{Q) and the Kesten-McKay support is smaller than aV~°', 
with a > and a positive constant a. Numerical simulations show that a ~ 2/3. 

• The cardinality of the set ^{Q) is bounded, and (|7(^)|) = C'(l). 

The derivation of the trace formula for Ramanujan graphs (2.14) used (2.7) 
which expresses the identity operator in terms of the orthonormal set of Chebyshev 
polynomials. Truncating the infinite sum in (2.7) aX, t — tmax, one obtains the identity 
operator in the finite space spanned by the first tmax Chebyshev polynomials. The 
resulting sum 

displays a peak centered ai x — x' of width The coarse-grained spectral density 

for fjL in the Kesten-McKay interval is defined as 

We return now to (2.5) with w = 1, and separate the sum on the spectrum to its 
Ramanujan and non-Ramanujan components: 

a = i (i - + ^ (<i - 1)-' (1 + (-1)') 

Upon multiplying both sides by ^^i^^^ o) '^^^ 2jl-i ^ ^^"^ performing the sum over 
t — 0, . . . , tmax we have to consider in particular the last term since it involves 
the exponentially increasing contributions from the non-Ramanujan spectral values. 
Consider ^ liQ)- It is bounded by < 2\/ri — 1(1 + aV^"") and therefore the 
contribution of the sum of the last terms in (2.19) is bounded by 

2|7(g)|exp(t^.,ay-"/^) 

V ^/ay-"/2 

Choosing t^ax < log V'^^'^^'^ ensures that the sum converges to zero as V increases. 
Summing the other terms in (2.19) results in the finite tmax analogue of the Kesten- 
McKay density, which converges exponentially quickly to the limit expression (the 
correction goes to zero as z^""''^, z < 1). The truncated oscillatory contribution arises 
from the truncated sum over the left hand side of equation (2.19). Truncating the t 
sum at tmax implies that the trace formula cannot resolve spectral intervals of order 
which are larger than the mean spectral spacing which is of order ^. In many 
applications this does not pose a severe problem. 
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2.2. The Magnetic Case (w^l) 

The magnetic spectrum differs from the corresponding non magnetic one in one 
important aspect, namely, there is no analogue to the phenomenon that the maximal 
eigenvalue //q is identically d. Rather, when averaged over the magnetic ensemble, it falls 
right on the boundary of the Kesten-McKay support. This can be shown by the following 
heuristic argument (A rigorous proof can be found in [43]). Consider the magnetic edge- 
connectivity matrix which excludes back-scatter: Y^^^ = — J. The ensemble mean 
of its maximal eigenvalue can be estimated by studying the behavior of ^tr (Y^^^Y ^ 
for large t, where (■ ■ denotes the average over the magnetic ensemble. The only non 
vanishing contributions to ^tr (Y^^^Y ^ come from self-tracing periodic nb- walks 
where each edge is traversed an equal number of times in both directions. The periodic 
walks which asymptotically dominate ^tr (Y^^^Y ^ consist of two cycles which share 
a single vertex, and each cycle is traversed twice in opposite directions. The common 
vertex enables the inversion of the traversal direction without back-scatter. On average, 
the number of such periodic walks is of order t{d—l)^ hence, {\r]o\)>s> — Vd — 1, where r^o 
is the largest eigenvalue of y(^) in absolute magnitude. Using the connection between 
the spectra of y(^) and A^^^ which is imphed by the Bartholdi identity, we deduce 
(|//o|)$ ^ 2Vd - 1. 

Starting from equation (A. 7), we continue by assuming that the graph is Ramanujan 
in the sense that all eigenvalues of A^^\ including the largest one satisfy < 2y/d — 1. 
The spectrum of Y^^^^ consists of 

^ {+ix{E-V), -lx{E-V), 

{Vd-1 e''t'\ Vd-1 e-*^^ A; = 1, • • • F)} 

where (j)k — arccos — , < 0^ < tt. (2.20) 

2v(c^- 1) 

The scaled traces yt are defined as 
1 tr(F(^))* 

yt = 77 — , (2.21 

and we end up with a trace formula, similar to equation (2.14). The smooth part is 
again the Kesten-McKay density. The oscillatory part is different, because the terms in 
yt contributed by t-periodic walks are now decorated by magnetic phases which are 
accumulated along the walks. We have shown previously that ^tr (f(^))* ^ ~ 

t{d — 1)^. Thus, {yt)<s> — >■ if the limits t ^ oo and V ^ oo are taken such that 
|f — > 0. Therefore the mean spectral density is again the Kesten-McKay measure, as in 
the non magnetic w — 1 case. 

3. Applications of the Trace Formula 



The w-trace formula offers a bridge between the spectral and the combinatorial aspects 
of graph theory. This connection can be exploited in order to compute combinatorial 
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quantities, as will be shown in the present section. 

Throughout this section we shall assume the large V limit of the spectral density 
(2.9), and neglect the term ^p'^°"'{fi;w): 

The explicit expressions for p^"^°°*^{p;w) and p°^'^{p]w) are given in (2.10). 

In the first section we have shown that ii Y^{w) can be considered as a generating 
function for the number N{t; g) of t periodic walks which scatter back exactly g times: 



N{t- g) 



'-1)9 Qg t^YHw 



g\ dw^ 



(3.1) 



=1 



Here, we shall use the w-trace formula to compute N{t,g = 1) explicitly and show 
how expressions for higher g can be obtained. To start, we take the first derivative of 
(2.9) with respect to w. The left hand side vanishes, since p{ii) does not depend on w. 
Moreover, it can be easily checked that '^^"'"""^ ^^'"'^ = 0. Thus, the first derivative of 
p°^'^{p; w) computed at w = 1 must vanish identically for any p. We shall show that the 
first derivative can be written in the form ^'-^K^^^^)' therefore the coefficients 
ai must vanish. This requirement provides a recurrence relation from which N{t, g — 1) 
can be computed for any t. 

From now on we shall denote ^ by ()'. Taking the derivative is quite tedious, since 
w appears in (2.10) both in the coefficients of the Chebyshev polynomials and in their 
argument. However, using elementary relations between the Chebyshev polynomials 
and their derivatives, and after some lengthy but straightforward computations, one 
gets: 

ai = {d- 2) Q(/ - 2)yi_2il) - ^(/ + 2)y^+2(l) - yz(l)) + 

{d-l)(^y[{l)-lyUi)-lyU^)) (3.2) 
We now recall that: 

VM - y ' . ^ (3.3) 

and we write (3.2) in terms of the try': 

U,.. ... 2)"-i-"(l)_(,_q(, + 2)"-^''<') 




(d-l)T (<i-l)" (rf-l)l J 

The requirement that = 0, results in an expression for (trY'')' in terms of {trY'^), 
where k < I. It is convenient to define: 

_ (try')'(i) 

„ ^ I (rf - 2)(i - 2)|^^ - (d - 2)(i + 2)^^ I (3.6) 
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and after some further computations, the following inhomogeneous recursion relation 
emerges: 

Efe=i92fc, if ^is even 



Pi+2 - Pi 



EI=o '?2fc+i, if ^is odd 
Substituting qi from(3.5), we get that: 



Pi+2 - Pi 



-2{d-2) 



-2{d-2) 



1/2 fry^fc(l) I trY\l) 



U/2J try2fc+l(l) ; fry'(l) 
k+h ~^ 2 



if / is even 



if / is odd 



((i-l)'=+2 " (d-l)2 

It can be verified by substitution, that the solution of the recursion relation is: 



l/2-l trY'^'^il) 

(d-l)'' ' 



if / is even 



Pi 



-Kd-2)j:^^Lf-rX^, if /is odd 



or, equivalently: 
(trF')'(l) = 
Finally, N{t; g 



Kd-2)Y:H-;UrY^\l){d-l) 



if / is even 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



-l{d-2)Y}i=t^trY^^+\l){d-l)'=^, if / is odd 
1) = — (try')'(l) is the number of periodic walks of length / with 




Figure 2. A periodic walk of length 16 composed of a cycle of length 10 and a tail of 
length 6 



one back-scatter. This result can be derived directly by the following argument. For 
simplicity, let us consider the case where / is even. A periodic walk of length / with one 
back-scatter must consist of a periodic walk of length 2k with no back-scatter, and a 
'tail' of length / — 2k. A tail is a periodic walk which goes over some path, and then 
comes back in the opposite direction. Clearly, it has one back-scatter, as demonstrated 
in figure (2). The number of periodic walks of length 2k with no back-scatter, is just 
trY'^^{\). The first edge in the tail can be chosen in c/ — 2 ways. Any other edge, until 
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we reach the end of the tail, can be chosen in d — 1 ways. There are such choices 

to make. The way back along the tail can be made in only one way. We now need to 
sum over k. Since both /, and the length of the tail, are even, the shortest periodic walk 
with no back-scatter is of length 4, so the sum must start at A; = 2, and it must end at 
k = I — 1, since the shortest tail is of length 2. 

We have thus seen that using the w-trace formula, we were able to extract 
combinatorial information about the graph. The computation above was the simplest 
case, but it is clear that N{t; g > 1) can be derived by further differentiations, and 
the use of known recursion relations involving the Chebyshev polynomials and their 
derivatives. 



4. Unitciry Evolution 

The trace formula (1.13) which was discussed in section (1) was originally derived by 
constructing the graph evolution operator U (jj) (1-14). The zeros of the secular function 
z{ii) = det(/ — U{ii)) provide the spectrum of the graph, and using standard methods, 
the trace formula followed. In view of the results derived in the previous sections, it is 
natural to expect that together with C/(/x), there exists a one parameter family of unitary 
operators, by which other secular equations can be written down, and corresponding 
trace formula can be derived. This is indeed the case, and to prove it, we go back to 
the Bartholdi identity for regular graphs (2.1). It is convenient to slightly modify the 
free parameters w, s and use a = s , j3 = ws, Bartholdi's formula reads now: 

det(/(2^) - {aB - (5 J)) = (1 - /J^)^"^ det((l + da(5 - (3^)1^^^ - aA) (4.1) 

for any complex a ^ 0, /9. 

Defining the IE x IE matrix U = aB — j3J and using the properties of B, J it is 
straightforward to show that 

' |/5P = 1 

\afd = a(3* + a*(3 



[/ is unitary ^ ^oth { \[\,^ „ ■ (4.2) 



The above implies that /3 = e**^ where is an arbitrary real parameter. If we impose 
further the relation 

a// = 1 + daP - 0^ (4.3) 

between a,j3 and the spectral parameter /i, Bartholdi's formula may be written as: 

det(/(2^) - U{n, 4))) = (1 - e2^^)^(/x - de'^)-^ det(/x/(^) - A) (4.4) 

where 

1 - e^''^ 
fi — de^ 

is a unitary evolution operator depending parametrically on the real and independent 
parameters and in the domain V = {/i E {—d,d), e M \ {ttZ}}. is restricted 
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away from zero and any integer multiples of tt since for these values a = and (4.4) 
does not provide a relationship between the spectrum of A and U{ii, 0) — — J. 
For ij,,4> the secular equation 

Za{ii) = det{iiI-A) = (l-e2^^)-^(At-de''^)^det(/(2^)-C/(/x») 

is the characteristic polynomial of A which is real on the real axis and its zeros coincide 

1^ 

V 



with the spectrum of A. The spectral density function, p{fj,) = ^ <^(a* ~ A*j)) which 



can be written as: 

pifj) = hm Jm log Za{ij> + ie) (4.5) 

VTT e->0+ djJ, 

will now be a sum of a "smooth" contribution from the phase of the (1— e^*'^)^^(/x— rfe*'^)^ 
term and a fluctuating contribution from the sum over periodic orbits on the graph, with 
amplitudes which arc determined by the forward and backward scattering defined by 
the evolution operator U (/i, 0): 

Pin) = - log ((1 - e'^^r^ip - de^^Y) 

1 d °° 1 

+ — lim Xm — V - tr iUip + ie)*) (4.6) 

V-K e-^0+ dpj^^t^^'' ^ ' 

Various real functions = (pip) can be defined, which yield various smooth and 
fluctuating parts. One interesting case is choosing (f) — const. For this choice the 
spectral density function takes the form: 

P^^^ ^-Vd (^)^ + l-2Mcos0 ^'-^^ 

1 d °° 1 

+ — lim Xm — V - tr (Vip + ie)*) (4.8) 
V-K 6-»0+ dp^^ t ^ ^ ' ' ^ ' 



where 

tr[/(^)* 



(//^ + — Ipd cos 0) 2 



N{g;t) aP ^ -</'+arctan 2l-i,cos4>+d 00^24, 1 (4.9) 



and 

. , J2 + d(d - 2) + u2 - 2pd cos + 2(d - 1) cos 20 

2sin0 

Speciflcally, the case of = — | corresponds to the evolution operator chosen in 
(1.13). Another interesting choice would be that for which the smooth part is the 
Kesten-McKay measure (1-12). This will be achieved for 0(//) which satisfles: 

d 2smo- + 2) -2i(l- l)//c()S(i>l 

V- 2^,2 0, \ = • (4.10) 

27r p^ + d^ — 2pdcos(p 
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where the L.H.S of this differential equation is the general expression for the smooth 
part of the spectral density function, normalized by the number of eigenvalues V, or 
equivalently 

V 2\2 n J 71 ^ 

:F- arccos ^^^^ - (-1 zp 1) = Nkm{i^) (4.11) 

TT \^J fj? -\- d^ -2 fjid COS (j){fji) j 

where the L.H.S is the general expression for the phase of (1 — e^*'^)~^(/i — de^'^Y 

divided by ttF, Nkm{ij) = ^ Pkm(a*) is the Kesten-McKay counting function and 

the =F refers to the sign of fism(f). The parameter A; e Z arises from the fact that 
(1 — e^*''^)~^(yU — de^"^)^ is defined up to an integer multiplicity of 2ti. Note that if 
'i^^KM^I^) solution to (4.11), then (p^j^j^j^jJi) + 27r is the k + ^'^^'^^^ solution, so 

in fact there are at most ^'^^^^^ distinct solutions to (4.10) and (4.11), where all the 
rest result by an addition of multiples of 27i. Figure (3) shows the numerical solution of 
(4.10). 

Assuming very large d and scaling the spectral parameter, u — ^^^rf' expand 
in powers of so that (4.10) takes the form: 

and solve to first order: 

AW)^.^ 1 . fjn T\ 2 fi 



™(/^) = Ck- wr, — rT^/i\/4(d- 1) - /i2 - - — - arcsin 



2(d-l)2-^ V ^ ^ d-i 2Vd^ 

where = | + + (4-13) 

To every distinct solution corresponds a set of continuous operators U^\^{iJ,) = 
U{n, 0^m(a'')) for which the smooth part is the Kesten-McKay measure (r = y). Thus, 
the fluctuating part must vanish after averaging over the ensemble, as V goes to infinity. 
Although a solution (j)^^^^ cross the value 0^m(a*o) = ttA; for /xq 7^ ±d, hence is not 
a valid solution under the construction above, a continuous set of unitary operators can 
still be defined using this solution by: 



U^'\lA = { ^-"""^'^^ (4.14) 



with the desired property. 



5. Discussion 



The novel feature which emerges from the present work is the fact that the same spectral 
density can be expressed by trace formulae which are based on different sets of periodic 
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Figure 3. The numerical solution of (4.10) for d = 4 and 



2k _ 
V ~ 



walks. Moreover, the form of the smooth part and the weight associated to periodic 
walks depend continuously on the parameters w and (f) in (2.10) and (4.6), respectively. 
This is most prominent in the case w = 1 where the contributing walks are non back- 
scattering. For any w ^ 1, back-scattering walks are contributing, so the weights and 
the phases associated to each periodic walk are w dependent. 

The choice w = 1 was beneficial for several reasons. First, it serves as a new and 
independent proof of the Kesten-McKay formula (1.12), where the tree approximation 
is not taken explicitly. Second, we were able to derive a leading order correction for 
large, but finite, V. 

The trace formulae we have derived, provide new insight into the question raised in 
the introduction: Can we distinguish amongst the periodic orbits on the graph distinct 
subsets, each responsible to a different feature in the spectral density? As we have 
mentioned, the smooth part (the Kesten-McKay measure) stems from periodic walks 
restricted to tree-like subgraphs. We can now identify the nb-periodic orbits as the 
subset which is responsible for the fluctuating part of the spectral density. It is therefore 
conjectured that for w 7^ 1, the back-scattered walks give opposite contributions to the 
smooth and oscillatory parts of p{fi). This sheds a new light on the puzzle regarding 
the re-summation which is required in (1.13). 

We have chosen to write the fluctuating part of the spectral density, (2.14) in a 
slightly peculiar way. Usually in trace formulae, the summands in the sum over periodic 
orbits, can be written as an amplitude times an exponent. The exponent plays the role 
of classical action. We can thus identify the quantity arccos ^ as an 'action' per 



unit step, of the quantum evolution on the graph. The amplitude is proportional to yt, 
which is a measure of the number of t-periodic nb-orbits. 

As we have shown, the w-trace formula can be used as a tool for turning spectral 
information about graphs into combinatorial information. This is one of the advantages 
of the w-trace formula over the one specialized for w = 1. 

In the last section, we have proven that a set of unitary evolution operators, which 
govern the quantum evolution on the graph, exist, and that they reproduce the density of 
states for the discrete Laplacian. This stands in contrast to the non-unitary "evolution 
operators" which were used in the preceding sections. As was shown, a specific choice 
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can be made, in order to reproduce the Kesten-McKay measure (or the semi-circle 
distribution, for large d). 

In the following paper in the series, we shall use the fluctuating part of (2.14) in 
order to investigate the spectral fluctuations of discrete Laplacians on regular graphs. 
We shall show that these fluctuations are given by the appropriate ensembles of Ran- 
dom Matrix Theory, and estabhsh some interesting connections between spectral and 
combinatorial graph theory. 

Appendix A. The Bcirtholdi Identity 

The basis for the work presented in this paper, is an identity proved recently by L. 
Bartholdi [20] which reveals a number of interesting connections between the spectral 
properties of the graph and its set of periodic orbits. This identity generalizes a previous 
result by Bass [21]. For general graphs, the Bartholdi identity is: 



det(/(2^) - s{B - wJ)) = (1 - iv\s^ f-^ dct(/(^) + ws\D - wl'^^^) - sA) . (A.l) 



identity matrices in dimensions 2E and V , respectively, and the matrices A. B, D and 
J were defined in the (1). For ri-rcgular graphs, the Bartholdi identity reads: 



det(/(2^) - s{B - wJ)) = (1 - wh^f-^ det(/(^)(l + w{d - w)s^) - sA) . (A.2) 



Its importance in the present context comes from the fact that it connects the spectrum 
of the adjacency matrix A with that of the matrices B and J, which can be used to 
count various types of cycles and walks on the corresponding graph. 

Bartholdi's original proof of this identity is based on combinatorial considerations. 
Proofs which are algebraic in nature were provided by various authors, and we quote 
here the version given by Mizuno and Sato [22] to make the present paper self contained, 
and to provide a base for the derivation of the trace formula for the discrete Laplacian. 

Appendix A.l. Proof of Bartholdi's identity 

The proof of Bartholdi's identity follows almost verbatim the proof presented in [22] . 
Define the two 2E x V rectangular matrices 



The parameters s and w arc arbitrary real or complex numbers, /(^^) and I^^^ are the 




(A.3) 



Denoting by X the transpose of X, one can easily prove that 

B^+^B(^) = B ; B(^)B^+^ = A 



(A.4) 
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Construct the two {2E + V) x {2E + V) square matrices 







-B(-) + wsB(+) 



j{2E) 



M 



ws 



Bi+) 



(1 - w;V)7(2£^) 



(A.5) 



Using the identities (A.4) one can compute the matrices LM and ML, and since their 
determinants are equal, one finally gets the desired identity (2.1). 



Appendix A. 2. Generalizations of the Bartholdi identity 

The Bartholdi identity can be easily generalized to three interesting cases: "Magnetic" 
graphs, multigraphs and weighted graphs. The proof always follows the same steps 
shown above, and for each case, one has to define the four matrices: A, B, D, J properly. 

• Magnetic regular graphs: 

We replace the adjacency matrix A by its "magnetic" analogue (1.10), and the edge 
connectivity matrix B by 

(sW),,e' = B,,,J^t>e+<P.'^ (A.6) 

The matrix J is not modified since 0e — ~0e- The "magnetic" Bartholdi identity 
now reads: 

det(/(2^)-s(5(^)-wJ)) = (l-wV)^-^det(/(^)(l+^^;(rf-^^;)s2)-sylW) .(A.7) 

The proof follows the same steps as above, after modifying -Bg^'' by multiplying them 
by e^5<^''=j and, by replacing the transpose operation (~) by hcrmitian conjugation. 
For non-regular magnetic graphs, the matrix D is the same as the non-magnetic 
case. 

• Multigraphs: 

The adjacency matrix is defined as explained in the introduction (section (1)). B 
is still a (0, 1) matrix where we must list all the edges, including parallel ones and 
loops. J does not change and in D we must count the degree of a vertex including 
parallel edges and counting loops as two edges. 

• Weighted graphs: 

The weighted adjacency matrix was defined in section (1). J is not changed, and 
we redefine B,D in the following way: 

(S(^))e,e'=Se,e'A/W^ (A.8) 

{D^'^% = 6,, (A.9) 

e:t{e)=i 

These are three canonical generalizations. Obviously, one can make further 
generalization by combining them (magnetic multigraph, for example). 
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